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BARYCENTERS OF MEASURES TRANSPORTED BY 
STOCHASTIC FLOWS 
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Universite de Poitiers and Nottingham Trent University 

We investigate the evolution of barycenters of masses transported 
by stochastic flows. The state spaces under consideration are smooth 
afline manifolds with certain convexity structure. Under suitable con- 
ditions on the flow and on the initial measure, the barycenter {Zt} 
is shown to be a semimartingale and is described by a stochastic 
differential equation. For the hyperbolic space the barycenter of two 
independent Brownian particles is a martingale and its conditional 
law converges to that of a Brownian motion on the limiting geodesic. 
On the other hand for a large family of discrete measures on suitable 
Cartan-Hadamard manifolds, the barycenter of the measure carried 
by an unstable Brownian flow converges to the Busemann barycenter 
of the limiting measure. 

1. Introduction. We consider the motion of a mass moving according 
to the law of a random flow. This can be used to model the motion of 
passive tracers in a fluid, for example, the spread of oil spilled in an ocean. 
Such motion can be assumed to obey a stochastic flow where particles at 
nearby points are correlated, for example, isotropic flows as investigated in 
[22, 27] or a general semimartingale flow as considered in [5]. The evolution 
of pollution clouds in the atmosphere or a gas of independent particles, 
on the other hand, can be described as blocks of masses moving according 
to the laws of independent stochastic flows. Here we propose to study the 
dynamics of masses transported by stochastic flows by investigating the 
motion of their centers of mass. As the medium in which the liquid travels is 
not necessarily homogeneous or flat, it makes sense to work on a nonlinear 
space, for example, on a manifold diffeomorphic to W n but with different 
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geometric structure. All measures considered here are normalized to have 
mass 1. 

The center of mass for a measure on IR n is the minimizer of the square 
distance function averaged with respect to the measure. This definition is 
used to define martingales on R n . The same minimizing problem can be 
considered on Riemannian manifolds. Although the distance in L 2 is the 
traditional object to minimize, we can also make sense of finding the min- 
imizer of the distance function in L . Busemann barycenters, for example, 
of measures on the boundaries of hyperbolic spaces are minimizers of Buse- 
mann distance functions in L . The Busemann barycenter was investigated 
for continuous measures for its application in analysis (see, e.g., [3]), while 
we are interested in both discrete and continuous measures. We shall re- 
late the center of a mass pushed by a stochastic flow, in the limit, to the 
Busemann barycenter of its limiting measure. 

From the point of view of minimizing assumptions on the state space, 
we note below that a center of mass can be defined using the concept of 
geodesies. For this we only need a linear connection V on a smooth mani- 
fold, and the pair (M, V) shall be called an affine manifold. Denote by 7„ 
the geodesic with initial velocity u S T X M . Denote by exp : TM — > M the 
exponential map exp^ u = 7 U (1), if it exists. 

A point x 6 M is an exponential barycenter of \x if it is a solution to 



In the sequel, by barycenter we always mean exponential barycenter. Note 
that (1.1) always makes sense locally and therefore for measures with suf- 
ficiently small support. However, to consider more general measures we are 
obliged to work on convex manifolds. An affine manifold (M, V) is said to 
be convex if for every pair of points x,y E M there exists a unique geodesic, 
defined using V, joining x and y and this geodesic depends smoothly on 
x and y. 

Main results. In Section 2 the existence and uniqueness of exponential 
barycenters of a probability measure with compact support are established 
for CSLCG manifolds. A CSLCG manifold is an affine manifold satisfying 
certain convexity conditions. Furthermore, a stochastic differential equation 
governing the motion of the barycenter is given. 

To study the large time behavior of the barycenter, we first investigate the 
toy model of the empirical measure \it = \ Ya=i °f n independent Brow- 
nian particles {X\) on an hyperbolic space HI. For n = 2 this is explicitly 
studied in Section 3 and for n > 3 this can be considered as a special case of 
the discussion in Section 4 and can be proved analogously. For two particles, 
when t gets large the barycenter gets close to a Brownian motion of variance 
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1/2 on the geodesic with asymptotic directions and X^, where X^ is 
the limit of X\. For more than two particles, the barycenter converges to a 
unique point in the manifold. This result seems to be surprising. However, it 
has a link to the elementary Steiner problem of finding a point which mini- 
mizes the sum of its distances to n given vertices: there is a unique minimizer 
if and only if the vertices do not all lie on the same line. If we replace the 
usual distance in this problem by the Busemann distance function and look 
for the point which minimizes the sum of the distances to the n given points 
on the boundary, we shall obtain the Busemann barycenter of the empirical 
measure of the limiting points. The unique point the barycenter of fit con- 
verges to is this Busemann barycenter. For n = 2, the two limiting points 
can be joined by a geodesic and every point on the geodesic is a minimizer 
for the Busemann distance function. 

In Section 4 the state space under consideration is a Cartan-Hadamard 
manifold with pinched negative curvature on which the Brownian motions 
satisfy a law of large numbers. In this case a Brownian flow Ft{y) converges 
almost surely for every starting point y to a process -Foo(y) on the visibility 
compactification of M . Furthermore, under suitable conditions, the limiting 
process separates the initial points and the uniform o{t) fluctuations for the 
Brownian motions vanish in the limit. Thus the L 2 minimizer of the distance 
function is the same as the L 1 minimizer, the Busemann barycenter, of the 
limiting measure on the boundary. More precisely, suppose fiQ has finite 
support and all the weights are smaller than 1/2. Then the barycenter of 
fit converges almost surely to a random variable 6 M , the Busemann 
barycenter of .Foo(^o)) characterized by 

lim f <p{Z oo ,F t {y)){0)d»{y) = 0, 

where ((p(z,y)(s),s > 0) is the unit speed geodesic connecting z and y start- 
ing at z and ip(z,y)(0) = -^\ s=0 tp(z,y)(s). 

Notation. Throughout the paper, a standard filtered probability space 
(£7, JF, (J r t)j>o,P) is fixed. If (Xt) is an M-valued continuous semimartin- 
gale and a is a C 1 section of T* M , we denote by J (a(X r ), 5X r ) and 
/o (a(X r ), (f 7 X r ), respectively, the Stratonovich and Ito integrals of a along X 
The second integral requires a connection V on M. If Xt(iv) takes its values 
in the domain of a local chart for any t G [S(u>), T(co)[ where S and T are 
stopping times, then 

j\a(X t ),5X t ) = j^a^dXl + \ ^ d -^L d {X\X% 

and 

j\a(X t ),d v X t ) = J s aiiXt^dXi + lr^X^diX^X^t), 



4 



M. ARNAUDON AND X.-M. LI 



where T l - k are the Christoffel symbols of V. The Ito integral can be defined 
via the Stratonovich integral as follows: let //o,t : Tx M — > Tx t M be the 
stochastic parallel transport along (Xt) and let Zt = Jq //q lSX r be the anti- 
development of (Xt), a Tx M- valued semimartingale. Then 

f\a(X a ),<Fx a ) = I (a(X s ) o // ,„dZ s ). 
Jo Jo 

A semimartingale (Xt) is called a V -martingale if and only if, for every 
smooth section a of T*M, ^(a(X s ),dy ' X s ) is a local martingale. Equiva- 
lently the stochastic anti-development of (Xt) by V is a local martingale. 
The reader may consult [7, 8, 9] and [12] for general reference. 

2. Exponential barycenters of measures transported by flows. 

2.1. Preliminaries. If the geodesies determined by two linear connections 
are the same, then the two connections differ by a tensor of type (1,2). In 
particular, given any connection, we can always subtract from it half of 
its torsion T to obtain a torsion- free connection with the same geodesies. 
Since convexity and barycenters are determined by geodesies, we assume 
that V is a torsion-free connection. We do not impose the condition that 
our connection be metric with respect to any Riemannian metric on M. 
However, for calculations we shall fix an auxiliary Riemannian metric on M 
with the corresponding Riemannian distance function denoted by p. In the 
sequel when the term " Riemannian manifold' is used it is considered as an 
affine manifold with respect to the Levi-Civita connection. A smooth local 
diffeomorphism between two manifolds M and (M, V) induces a connection 
V on M. An affine map (j> between two affine manifolds (M, V) and (M, V) 
is a smooth map such that for all smooth functions / : M — ► R, Vd(f o <fi) = 
4>*(Vdf). An affine map on M preserves geodesies. We shall frequently use 
the following: 

1. A point x in M is a barycenter of fi if and only if 

(2.1) H(x) = H„(x) = [ 7(x,y)(0)M<%)=0 

Jm 

for geodesies ( r y(x,y)(s) :0 < s < 1) connecting x and y. 

2. If z is the barycenter of \i and G:M — > N an affine map between two 
convex affine manifolds, then G(z) is the barycenter of the push forward 
measure G(fi) by G. 

3. For a smooth Riemannian manifold let fi be a probability measure on M 
such that f(x) = J M dist 2 (x, y) fi(dy) is finite for some (hence for all) x 
in M. Then the exponential barycenters of /i are the critical points of 
the function / ([10], Proposition 3). Moreover, if x is an exponential 



BARYCENTERS OF MEASURES 



5 



barycenter of p and h is a convex function on M, then we have Jensen's- 
type inequality ([10], Proposition 2) 

(2.2) h(x)< f h(y)n(dy). 

JM 

A convex manifold M is diffeomorphic to an open set of R m , m = dim(M) , 
and there exists a neighborhood U of the null section in TM such that 

(a;, u) G U h-> (x, exp x u) G M x M 

is a diffeomorphism. Standard convex complete Riemannian manifolds in- 
clude Cartan-Hadamard manifolds (complete, simply connected manifolds 
with sectional curvature less than or equal to 0). Examples of incomplete 
convex manifolds are given by: small geodesic balls in Riemannian manifolds 
or small balls centered at the origin in an exponential chart in an affine man- 
ifold. [Note that every point x in a Riemannian manifold M has a convex 
geodesic ball B r (x); in the case of a sphere of radius r, a maximal convex 
set is given by an open half sphere, and its convexity radius is nr/2.] 

Recall that a function (f> : M — > R on an affine manifold (M, V) is convex 
if (f> o 7 is a convex function for all geodesies 7, that is, for all geodesies 
(7(x,y)(t),0 < t < 1) connecting two different points x and y, 

(2.3) <f>(rf(x,y){t)) < tc/>(x) + (1 - t)<j)(y). 

For a C 2 function </>, it is convex if and only if ^4t ( f ) ('Jt) — 0, where § 
denotes covariant differentiation with respect to t. 

For the uniqueness of exponential barycenter we often use a convex func- 
tion on M x M which separates points. The definition below is partly bor- 
rowed from [9]. 

Definition 2.1. Consider an affine manifold (M, V) and M x M with 
product connection. A convex function <f>:M x M — > R + vanishing exactly 
on the diagonal A of M x M is called a separating function on M. Let 
p £ 2N. A manifold which carries a smooth separating function (j) such that 

(2.4) eff < <p < CfP 

for some constants < c < C and some Riemannian distance function p is 
called a manifold with p- convex geometry. 

For Cartan-Hadamard manifolds, the Riemannian distance function is a 
separating function. Note that the square of the distance function is also a 
separating function and it is smooth. In more general Riemannian manifolds, 
sufficiently small geodesic balls have 2-convex geometry. For instance, any 
geodesic ball strictly smaller than an open hemisphere has p-convex geom- 
etry for some p depending on the radius, as proved by Kendall in [15], and 
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any point in an affine manifold has a convex neighborhood with 2-convex 
geometry. However, an affine manifold carrying a separating function is not 
necessarily convex. For example, M m \ {0} is not convex even though the 
distance function is a separating function. It is more difficult to find convex 
manifolds which do not carry separating functions; see [16]. 

Let (p be a C 2 convex function with {tp < 0} a relatively compact set. 
Then any probability measure [i on M with support included in {ip < 0} 
has an exponential barycenter, see [10, 14]. Also if M is an affine manifold 
with a bounded separating function, then any probability measure on M 
has at most one exponential barycenter [10, 14]. 

As a consequence: let M be a convex manifold with convex geometry, and 
let 4> be a C 2 separating function with 0~ 1 ([O, a[) relatively compact in M 
for some fixed x E M and some a > where <f> x = (f)(x, ■) (which is not the 
case in general). Then any measure supported on the set </>~ 1 ([0,a[) has a 
unique exponential barycenter. 

Observe that an open convex subset of a convex manifold with p-convex 
geometry is a manifold with p-convex geometry. 

2.2. Existence and uniqueness. 

Definition 2.2. A convex affine manifold (M, V) is said to be CSLCG 
(convex, with semilocal convex geometry) if every compact subset K of M 
has a relatively compact convex neighborhood Uk which has p-convex ge- 
ometry for some p 6 2N depending on K. 

Equivalently, a convex affine manifold (M, V) is CSLCG if there exists an 
increasing sequence (U n ) n >i of relatively compact open convex subsets of 
M such that M = Un>i U n , and for every n > 1, U n has p-convex geometry 
for some p £ 2N depending on n. 

This definition is motivated by Propositions 2.4 and 2.7 which state that 
any probability measure in a CSLCG manifold, of compact support, has a 
unique barycenter and that in a certain sense, the exponential barycenter of 
fi is differentiable as a function of \i. 

Note p increases with K. Examples of CSLCG manifolds are open hemi- 
spheres endowed with the Levi-Civita connection (which do not have p- 
convex geometry for any p £ 2N; see [15]). Examples of geodesically com- 
plete CSLCG Riemannian manifolds should be given by manifolds with a 
pole under curvature conditions to be determined. We conjecture that if 
(M, V) is a CSLCG manifold, then (TM, V c ) is a CSLCG manifold, where 
V c is the complete lift of V. Note that on every relatively compact convex 
subset of TM one can construct a continuous separating function: by [17], 
uniqueness of martingales with prescribed terminal values implies the ex- 
istence of such a separating function, and by [1], uniqueness holds. So the 
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conjecture concerns the smoothness of the separating function, and the fact 
that it satisfies (2.4). 

One can find in [16] an example of convex manifold which has not semilo- 
cal convex geometry. In the manifold constructed there, there exists a prob- 
ability measure carried by three points, which possesses four exponential 
barycenters. Clearly there is no convex neighborhood of these seven points 
in which we can define a separating function. 

Lemma 2.3. Let (M, V) be a CSLCG manifold. Every compact convex 
subset K of M has a convex neighborhood U with a nonnegative C 1 convex 
function cf>K such that (/^({O}) = K . 

Proof. Let U' be an open convex relatively compact neighborhood of 
K with p-convex geometry for some p G 2N. Let ^ be a smooth separat- 
ing function on U' satisfying cff < (J) <C ' fP where < c < C, and p is the 
Euclidean distance induced by a global chart, for example, an exponential 
chart. Define for x £U' 

(2.5) 4> K {x)=mf{4>(y,x),y£K}. 
Clearly ^({0}) =K, since K is compact. 

1. First we show that 4>k is convex on U' . Take {x\,X2) G U' x U' . Let 
p(xi) be points in the compact set K achieving the minimum: <j)K{, x i) = 
(j)(p(xi),Xi), i = 1,2. Since (j) is convex, for every t G [0, 1], 

4>(-/(p(x 1 ),p(x2)){t),-f(x 1 ,X2){t)) < (1 - t)(f)(p(xi),Xi) +t(/)(p(x2),X2). 

On the other hand, by the definition of 4>k, we have since j(p(x\),p(x2))(t) G 
K 

(t>K{l{xi,x 2 ){t)) < 4>(<y(p(x 1 ),p(x2)){t),"f{xi,X2){t)). 
Putting these equations together we obtain 

<f>K(Tf(xi,x 2 )(t)) < (1 -t)4> K (xi) +t(p K (x 2 ), 

which proves the convexity of (j>K- 

2. Next we show that (f>K is C 1 on some convex neighborhood U of K 
included in U' . We first prove that there exists a neighborhood U C V of K 
on which there is a unique point p{x) such that <I>k{%) = <f>{p(x),x). 

Suppose for x G U' and yo,yi G K, 4>k{x) = <f)(yo,x) = <f)(yi,x). Set y{t) = 
7(2/0 1 2/1 )(*)• Then necessarily, for every t G [0, 1], <fi(y(t),x) = 4>k(x)- 

By the Hadamard lemma (see, e.g., [4], Corollaire 3.1.9 and Exercice 3.1.10) 
we can write in the global chart 
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m p+1 

+ b h-i P+1 (y,x)~[[(y ij -x% 

h,...,i p+ i=l j=l 

where a^...^ and are smooth functions. 

Write f(t) = cj){y{t), x). It is a constant function. On the other hand, we 
can differentiate p times the function / with the expression of <f> in the chart. 
Using 

m 

y\t) = - ]T r jk (y(t))y J (t)y k (t), 
j,k=i 

where T l - k are the Christoffel symbols of V in the chart, we see that /W is 
of the following format: 

m p 

/«(*) = £ a H ... ip (x)l[y^(t)+g(y(t),y(t),x), 
il,...,i p =l j=l 

where g is a smooth function. Since U' is a relatively compact convex open 
neighborhood of K, we see by the Hadamard lemma and explicit calculation 
that 

\g(y,z,x)\ <C"||y-x||||z|| p 
for some positive constant C . Now since cff < 4>, we have for every z £ M. m 

m p 

E n ^' > 

h,...,i p =l j=l 

We choose U of the form U = {x' 6 [/', c^(x') < e}, with e G ]0, c p+1 /(C') p [, 
such that C7 is convex and relatively compact in U' . Consequently for x £ U 
and y £ K, \\y — x\\ p < (j)(x, y)/c < e/c and 

f^(t)>c\\y(t)\\ P -C'\\y(t)-x\\.\\y(t)\\ P 

>(c-(e/c) l ^C')\\y(t)\\ P Vi€[0,l], 
where c - {e/c) l l p C > 0. But / is constant, so we have y(t) = and conse- 
quently 2/o = 2/i- 

3. For x £U, we let be the point in K such that 4>k{x) = 4>(p(x),x). 
We prove that p is continuous on U. If not, let (x n ) ne N be a convergent se- 
quence in U with limit x such that p(x n ) does not converge to p(x). Since K 
is compact, by choosing a subsequence if necessary, we can assume — 
y £ K \ {p(x)}. By the continuity of 0, (p K (x n ) = 4>(p{x n ),x n ) -> 4>(y,x). 
On the other hand, since is convex on the open set ?7, it is continuous 
(see, e.g., [11], Proposition 1) and so lim^(a; n ) = 4>k{x). Consequently, 
4>(y,x) = 4>k(x) and by the uniqueness given in (2) we obtain y = p(x), a 
contradiction. So the map p is continuous. 
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4. We are left to prove that <ftx is C 1 . Denote by d g 4>K the Gateaux- 
differential of eft. Recall that for every x G U and v £ T X M, 

, , ( w n v ^(exptu) -(j> K {x) 
d g cft K {x){v) = hm 

and that dgcftxix) is convex on T X M (again by [11], Proposition 1). Let x £U 
and v £ T X M. If t > is sufficiently small, we have 

4>K (exp to) — (x) < 4>(p(x) , exp to) — <ftip{x) , x) 

which yields 

d g 4> K {x)(v) < d(ft p{x) (x)(v), 

(ft y denoting the map cft(y,-). Since d(ft p ^{x) is linear, this inequality im- 
plies d g (j)K{x) = d(ft p ( x -)(x) on T X M. Otherwise we would have d g cftjc(x)(v) < 
d(ft p ( x )(x)(v) for some v £ T X M, which would give by convexity of d g (j)K(x) 

d g <ft K {x){-v) > -dgcft K (x){v) > -d(ft p ( x )(x)(v) = d<p p ( x )(x)(-v), 

a contradiction. So dgcftxix) = d(ft p ^{x) on T X M . The differentiability of (ftx 
then easily comes from the inequalities, for all v £ T X M sufficiently close to 
0, 

< <ft K (expv) - (ft>K(x) - d(ft p ( x )(x)(v) 
< 4> P (x)(expv) - (ft p{x )(x) - d(ft p{x) (x){v). 
That (ftx is C 1 comes from the continuity of p. □ 

Proposition 2.4. Let M be a CSLCG manifold. Every probability mea- 
sure fx on M with compact support has a unique exponential barycenter. 

Proof. The proof is a slight modification of that of Proposition 5 
in [10]. For uniqueness let x and x' be two exponential barycenters of fi. 
Set K = supp(/i) U {x, x'}. Let eft be a separating function defined on a con- 
vex relatively compact neighborhood of K. Writing v the measure \x pushed 
forward by y \— > (y,y), (x,x') is an exponential barycenter of u, so 

ff(x,x')<-(ft(x,x')<- / (ft(y,y')du(y,y) = - / (ft(y, y) dpi(y) = 
c c J c J 

giving x = x' . 

For the existence let K be a convex compact subset of M containing the 
support of \x (we know that this exists in a CSLCG manifold). By Lemma 2.3 
there exists a C 1 convex nonnegative function (ftx defined on a relatively 
compact open neighborhood U of K, such that (ftJ^dO}) = K. Let e > 
satisfy (ftj}([0,e[) C U. We apply [10], Proposition 5, to the function (ftx — 
e/2 to see that /j, has an exponential barycenter in ^ 1 ([0,e/2[) (note [10], 
Proposition 5 is still valid with a C 1 convex function instead of a C 2 convex 
function) . □ 
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2.3. Mass moved by a smooth vector field. Let (M, V) be a smooth 
convex affine manifold with V torsion free. As before, for x, y in M let 
( , y(x,y)(s),0 < s < 1) be the geodesic, with respect to V, with end points 
x and y. The geodesic 7(x,y) shall be abbreviated as 7 where there is no 
risk of confusion. Set 

J(u,v)(s) =T(j(-,-)(s))(u,v), u G T X M, veT y M, sGK, 

where T-y(s)(u, v) denotes the derivative of j(s) in the direction of (u, v ). The 
map ( J(u, v)(s), < s < 1) is the Jacobi field satisfying boundary condition 
J(u,f)(0) = u and J(u,v)(l) = v. Write J(u,v)(s) = (V d / ds J(u,v))(s). For 
every x,y G M, we define the linear map 

(2.6) t/; (Xjy y.T x M^T x M, u » j(u, y )(0), 

where 0^ is the zero tangent vector in T y M . 

On TM one can define a connection V c , called the complete lift of V 
(see, e.g., [26]). It is a torsion- free connection since V is, so it is charac- 
terized by its geodesies, which are the Jacobi fields of the connection V. 
The canonical projection ir : (TM, V c ) — > (M, V) is affine. We easily see that 
(TM, V c ) is convex if (M, V) is and the only geodesic joining u and v in 
TM is ( J(u, v)(s), < s < 1). Given a C 1 vector field A, consider the map 
A : M — >■ TM. It induces a measure ^4(/i) on TA/. A point u G TM is a 
barycenter of a measure A(^) if and only if ir(v) is a barycenter of fi and 

(2.7) / j(v,A(x))(0)dfi(x) = 0. 

JM 

A measure // is differentiable along a vector field A if its pushed forward 
measure 5t(/u) is weakly differentiable in t, where for x G M, Sf(x) is the 
integral curve of A starting from x. The following lemma asserts that the 
map fi — ► 6(/i), 6(/i) the barycenter of /i, is a differentiable map in the above 
sense. 

Proposition 2.5. Suppose (M, V) is a CSLCG manifold. For i = l,...,n, 
let fi l be a probability measure on M with compact support K , let A 1 be a 
smooth C 1 vector field and let (Sl(x),0 <t<r) be its integral curve starting 
from x. For Pi> with J27=iPi = 1> se ^ = Pi (/•**)■ U z (t) * s 
exponential barycenter of fit : 

(a) then z{— ) is differentiable att = and i(0) is the exponential barycen- 
ter of Yl^iPiA 1 ^ 1 ) with respect to the complete lift V c of V; 

(b) there is Ck > depending only on K and the arbitrary metric such 
that 

p(0)||<C*r sup ||^(x)||. 

l<i<n 
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Proof. Let M' be a convex relatively compact open neighborhood of 
K with p-convex geometry (p 6 2N). Let be a smooth separating function 
on M' satisfying cff < (f> < Cff where p is a Riemannian distance func- 
tion, < c < C. For i < r the barycenter z(t) exists and by the definition of 
barycenters, (z(0),z(t)) is the exponential barycenter of YA=iPi(^hi ^Dit^) 
with respect to the product connection. 

Apply (2.2) to see <f>(z(0), z(t)) < Tl}=iPi Ik HS l (x), S l t (x)) dp\x), which 
implies 



p(z(0),z(t))< 



{21 



< 



-</>(z(0),z(t)) 
c 

C\ 1 / p ( n 



i/p 



i/p 



ff(S^x),Si( X ))dp\x) 



l - 



) t sup sup || A 1 (SI (x) 

J ie{l,...,n} x£K 



0<r<t 



where || • || is the Riemannian metric corresponding to p. Since all Rieman- 
nian metrics on M 1 are comparable, 



(2.9) 



1 i f t\ 

-exp z{Q) z(t) 



<^p(z(0),z(t)) 



<c o[ - 



C\ 1/p ■ ■ 

) sup sup 1 1 A 1 (S* (x)) ||, 

/ ie{l,...,n} x&K 



0<r<t 



some constant cq and so the family { | expj^ z(t), t < r} is bounded and 
thus has a limit point which shall be denoted by u. Let tk be a sequence 
going to such that lim^oo ^exp~^ z(tk) = u. Each zt k satisfies 



n , 

0*(t*)=Ett / j(z(tk),St k (x))(0)dS 
i=i JM 



and the covariant derivative of 7(z(i), Sl(x))(0) along the subsequence 
satisfies 

D 



dt k 



i(z(t k ),SUx))(0) = J(u,A i (x))(0). 



Integrate the identity over M and use the dominated convergence theorem 
to see 



n . 
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which together with 

n „ 

0,(o) =$> / j(z(0), x)(0)d^(x) 

fr[ Jm 

shows that u is an exponential barycenter of Y^iPi-^ifJ?) with respect to 
V c . 

Let 7T : TM — ► M be the canonical projection. To prove the uniqueness 
consider the map 

T® p (j>:TM' ->R+, lim — ^(tt(«), exp to), 

which is V c -convex and T® p <j){v) < C\\v\\ p . Furthermore the map 

(v, v') i— > T® p cj)(v' — v) 

is convex on £7 := {(■!;, i/) £ TM' x TM' , ir(v) = tt(V)} for the connection 
V c (g) V c (note £ is a totally geodesic sub manifold), and vanishes exactly on 
the diagonal of TM' x TM'. See [1] for the proof and details. 

Let u' be another exponential barycenter of A(p) with respect to V c . 
Since fx has a unique barycenter and it is affine, tt(u') = n(u). So it makes 
sense to estimate T®P(j)(u' - u). Set v = Y£=iPi{A\ A 1 )^). Then (u,u') is 
an exponential barycenter for u, so 

T® p (f>(u-u')< I T® p c/)(v-v')dv(v,v') 

JE 

n ,. 

= Ypi T® p <t>{A i {x)-A\x))dn i {x) = {). 

U. Jm 

This implies that u = v! , which in turn implies that the family | exp~^ z(t) 
has only one limit point as t goes to and so it converges to u. Thus 
1 1 — > z(t) is differentiable at t = with derivative the exponential barycenter 
of E?=iPi^- i (A ji ) with respect to V c . Finally, using (2.8) and (2.9), letting t 
go to 0, invoking the compactness of K and the continuity of S l , we get 

C\ l i p 

□ 



|z(0)|| < C K sup ||^(x)|| where C K = c 

ie{i,...,n} V c 

x<=K 



The map defined in (2.6) is not invertible in general, as can be seen 
with two points x and y at a distance 7r/2 in a hemisphere. However, we 
have the following: 



Lemma 2.6. Let M be a CSLCG manifold, let K c M 6e a compact 
set and let p be a Riemannian metric on M . Then there exists a constant 
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C' K > such that for any probability measure p in M with support included 
in K, letting z be the exponential barycenter of p,, the linear map on T Z M : 



ip {Z:X) (u) dp{x) 



M 

is invertible with inverse bounded by C' K (with respect to the metric p). 

Proof. Since M is convex, the map v i-» J(O a ,u)(0) from T X M to T Z M 
is injective and so is surjective. Its inverse — A(x, •) is well denned. For each 
u G T Z M fixed, A(-,u) is the vector field in M characterized by J(O z ,A(x, u))(0) 
—u. Let b(u) G T Z M be the exponential barycenter of A(u)(fi) (by Proposi- 
tion 2.5, it is well defined). By (2.7) 

= / j(b(u),A(x,u))(0)p(dx) 
Jm 

j(6(u),0)(0)/i(ds)+ / J(0,A(x,u))(0)p(dx). 
m Jm 

By the definition ij)r z ^ x \ (v) = J(v,0 x )(0) for v G T Z M, 



i>(z,x)(-)dn(x))(b(u)) = - / J(O z ,A(x,u))(0)dp(x) = u. 

M J JM 

We have sup xgA - \\A(x, u)\\ < C\\u\\ for some constant C since K is compact. 
Apply Proposition 2.5 to see 

(2.10) < C K sup \\A{x,u)\\ < C K C\\u\\. 

Thus f M ijj( Zi x)(') dp(x) is invertible with inverse satisfying 

(J il> M (-)dv(x)) (u) < ||6(u)|| <C^||u|| for C k = CC K . n 

As a direct consequence of Lemma 2.6 and its proof, we make more precise 
the result of Proposition 2.5 and give an expression for z(t). 

Proposition 2.7. Assume (M, V) is a CSLCG manifold with the con- 
nection V torsion free. For n > 1 and % = 1, . . . , n, let pt be positive numbers 
satisfying J27=i V% = ^- an d ^ A* 1 be probability measures on M with compact 
support. Let 

5?(-):[0,oo[x supp(^)^M 

be a map of class C k, ° for some k>l [i.e., C k in t with derivatives up to or- 
der k with respect to t continuous in (t,x)]. Then the exponential barycenter 
z{t) of pt := J2i'=iPi^t(l J ' 1 ') exists and is unique. Furthermore: 
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(i) t\-*z{t) is C k ; 

(ii) the bary center z(t) is characterized by the identity 

n . n . 

i=i jm i=i jm 

where Z is the zero vector in T Z M ; 

(iii) z(t) solves the ordinary differential equation 

i(i) = -|> j M {^Pjjj (zm{x) ) d ^^ 

(2 ' 11} xj(Oz(t),&t(v))(0)dlS(v), 
where ip is defined by (2.6). 

Taking (i % to be Dirac distributions, we get from Proposition 2.7 the fol- 
lowing: 

Proposition 2.8. If (M, V) is a CSLCG manifold with the connection 
V torsion free letn>\ andp= (pi, . . . ,p n ) be an n-tuple of positive numbers 
satisfying J27=iPi = 1- Then the map G p defined by 

G p : M n M, (x 1 ,...,x n )^b(f2Pi S x\ 
is smooth with derivative, 

n I n \ — 1 

TG p {ui,...,u n ) = ~YlPi\ j(0G(x u ...,x n ),Ui)(0), 
i=l \ j=l ) 

where (u 1: . . . ,u n ) £T Xl M x ••• x T Xn M . 

Proof. Just note if S l t (xi) = exp(W), then i(0) = TG p (u l , . . .,u n ). □ 

We want to generalize Proposition 2.7 to the case where for every x the 
flow is a semimartingale. We begin with a simple case, which is an immediate 
corollary of Proposition 2.8 since Z% as a smooth function of semimartingales 
is a semimartingale and stochastic calculus applies. 

Proposition 2.9. Let (M, V) be a CSLCG manifold. Let n > 1, p = 
(pi, . . . ,p n ) be an n-tuple of positive numbers satisfying Ya=iPi = 1> 
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be n points of M and /io = J2i=iPAi- F° r every i, let X\ be a semimartin- 
gale started at Xj. Then the exponential barycenter Z% of Y^iPi^x' ^ s a 
semimartingale and satisfies 



When /io is carried by two points we do not need to assume that the man- 
ifold has locally convex geometry: convexity is sufficient. If /io = \ {&x + <5 yo ) , 
Xt = Ft(xo), Yt = Ft(yo) for a stochastic flow Ft, then Zt is a semimartingale 
and 



(2.12) 5Z t = -^ {ZuXt) + ^ ZuYt) )-\j{Q Zt ,8Xt){Q) + J(0 Zt ,5Y t )(0)). 



2.4. Mass pushed by a random flow. Next we consider the case where Ft 
is a local semimartingale flow of homeomorphisms (which we abbreviate as 
semimartingale flow). Following Kunita ([18], Section 4.7) with small modifi- 
cations, see also [20, 21], let a(t,x,y,uj) be a predictable process with values 
in the tensor product T X M ®T y M , let b(t,x,uj) be a predictable process 
with values in T X M and let At be a continuous adapted real-valued increas- 
ing process. We say that Ft(x) is a semimartingale flow with characteristic 
(a(t,x,y,Lo),b(t,x,Lo),A t ) if 



dF t (x) ®dF t {y) = a(t, F t {x), F t {y)) dA t Vx,y G M, 
and the drift of Ft(x) is b(t, Ft(x)) dAt, that is, for every 1-form a S T*M, 



is a local martingale. 

Proposition 2.10. Let (M, V) be a CSLCG manifold with V torsion 
free and let /io be a probability measure with compact support. Assume Ft{x) 
is a semimartingale flow with characteristic (a(t,x,y,uj),b(t,x,uj), At), such 
that a(t, — , —,cu) and b(t, —,to) are C 1 with derivatives a.s. locally uniformly 
bounded in time and space. Then the exponential barycenter Zt of fit exists 
and is unique. Furthermore: 

(i) the process Z t is a semimartingale on [0,r[, where t is an almost 
surely positive random explosion time; 




If X\ = Ft{x l ) where Ft(x) is a semimartingale flow, then 
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(ii) the barycenter Z± is characterized by the identity 

[ j(5Z u Ft{x) )(0)dp (x) = - [ j(0 Zt ,5F t (x))(0)dp (x); 
Jm Jm 

(iii) Zf solves the stochastic differential equation 

(2.13) SZt = ~ J (J AzuF t {x))d»o(x)) j(0 Zt ,8F t (y))(0)d»o(y). 

Remark 2.11. The assumptions on the semimartingale flow can be 
weakened. We can assume, instead, that (a(t, x,y,w),b(t, x, u), At) belongs 
to the class -B ' 1 in the sense of [18], adapted to a manifold. 

Proof of Proposition 2.10. Fix an arbitrary Riemannian distance p 
on M. For simplicity we choose a and b bounded (by [18], remark on page 85, 
this amounts to changing At). Using a change of time, we may assume that 
At = t. With the assumptions on the local characteristics of the flow Ft(x), 
we can choose a version of Ft(x) jointly continuous in (t,x). Let (U n ) n >i 
be an increasing sequence of relatively compact convex open subsets of M 
containing supp(^o) and such that \J n U n = M . Define 

r n = inf{i > 0,F t (x) U n for some x G supp(/io)}. 

Then since the map (t,x) t— > Ft(x) is almost surely continuous and po has 
compact support, r n is an increasing sequence of stopping times converging 
to the explosion time r, which is almost surely positive. 

The map (£, y) i— ► Ft(y)(uj) is almost surely uniformly continuous on [0, R] x 
supp(//o) for every R > smaller than r. As a consequence, the barycenter 
Z t of pt is defined on [0,r[, and is continuous and adapted. 

We shall show that Zt is a semimartingale and is given precisely by the 
equations in parts (ii) and (iii) up to each time r n . Letting n tend to infinity 
will prove that the conclusion of the proposition holds. 

For the local result we begin with the following lemma which can be 
considered as an extension of Lemma 2.6. 

Lemma 2.12. Let U be a relatively compact convex open subset of M 
containing supp(po) and 

t\j = inf{t > 0,Ft(x) £ U for some x E supp(//o)}- 

There exists e > such that if t<T\j and p(z, Zt(u:)) < e, then 

/ 4>(z,F t (x)) dp (x) 

Jm 

is invertible with inverse bounded by a constant depending only on U . 
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Proof. Consider the set (U)' of probabilities on U, endowed with the 
weak topology. It is a metrizable space. Since U is compact, the map 

/ ip (zx) dp(x) 
Jm 

is continuous on U x (U)' for the product topology. Furthermore, the map 
from (£/)' to M sending p to its exponential barycenter z(p) is continuous 
by Skorokhod theorem: if p k converges to p, there is a sequence of random 
variables X k and a random variable X with values in U such that the law 
of Xk is /ifc, the law of X is and converges almost surely to X. Since 
U is compact, p(Xk,X) converges to in LP for every p 6 2N. On the other 
hand, (z(pk),z(fj,)) is the exponential barycenter of the law of (X k ,X) and 
U is p-convex for p sufficiently large. So 

p(z(p k ),z(p))<C\\p(X k ,X)\\ p 

which implies that z(p k ) converges to z{p). This proves that p \— > z(p) is 
continuous on (U)' . 

Next recall that for every p G (U)', J M dp(x) is invertible. The 

image of the compact set (U)' by the continuous map p\— > $m^{z(ij)>b) dp{x) 
(w.r.t. the weak topology) is compact. This implies the existence of a neigh- 
borhood of the graph of p i— > z{p) of the form {(p, z) € ((U)' x {/), p(z, z(p)) < 
e} for some e > 0, on which f M i/)r z x \dfj,(x) is invertible with inverse in a 
compact set depending only on U. This proves the lemma. □ 

We continue with the proof of the proposition. For relatively compact set 
U n , let e n be the constant defined by Lemma 2.12 and set 

T x n := r n A inf jt > 0, p(Z , Z t )>^. 

For every z £ U n with p(Zq,z) < e n /4 we define a T z M-valued stochastic 
process: 

G(z,t) = -£j (f 4(z,F,(x))d(M)(x)) j{0z,SF s {y))(0)dp {y), 

0<t<Tl 

Clearly G(z,t) is a semimartingale for each z with bounded local character- 
istics depending smoothly on the spatial parameter z (see, e.g., [27], Theo- 
rem (3.3)). So one can solve the equation 

(2.14) 5Z' t = G(Z' t ,5t), Z' = Z . 

Its solution Z[ exists and is a semimartingale up to time 

T' n := T x n A inf jt > 0, p(Z' t , Z t ) > ^ j . 
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We only need to show that Z[ = Z t , the exponential barycenter of Ft(fXo), 
that is, 

(2.15) / >y(Zl,F t {x))(0)di«i(x) = O zl . 
Jm ' 

The equality holds for t = 0. The process ■j(Z' t , Ft(x))(0) is a semimartin- 
gale with continuous spatial parameter and local characteristics uniformly 
bounded on (x,t) £ supp(//o) x [0, TJ^[ and so using [27], Theorem (3.3), we 
see that integration over M commutes with covariant Stratonovich differen- 
tiation Dt in t: 

(2.16) D t f j(Z' t ,F t (x))(0)d fM) (x)= [ D t j(Z' t ,F t (x))(0)dMx)- 

Jm Jm 

More precisely, (2.16) is equivalent to: writing Vt(x) = j(Zt, F t (x))(0), 
J* (p* (a), 6^ J m V s (x)dno(x)j^ = J m ( (p* (a) ,6V s {x)fj dfi (x) , 

for every < t < T' n and every section a of T*M. Here p:TTM — > TM is 
the map induced by the connection which to an element of TTM associates 
its vertical part. 

Now since the connection is torsion free and since 7 is smooth, -C^ = 
where Sr- denotes covariant differentiation in s. Consequently, we have on 

[ Da(Z' t ,F t (x))(0)dfi (x) 
Jm 

= f ^-6a(Z' t ,F t (x))(0)dfi (x) 
Jm ds 

= [ j(8Z' t ,8F t {x)){ti)dn Q {x) 
Jm 

= f j(6Zl,0 Ft(x) )(0) + j(0 z ,,6F t (x))(0)d» o (x) 
Jm ' 

= ( Jm ^( z 'f Ft ( x ^^ ^( x )) ( SZ t) + j M j (°z t ' . SF t (x))(0) d(i (x) . 
Plugging the expression (2.14) for bZ[ in the last formula yields 

(2.17) D t f j(Z' t ,Ft(x))(0)dno(x)=O z , 

Jm * 

Together with Z' Q = Zq we see that for t £ [0, T^[, 

(2.18) / i(Zl,F t (x))(0)diMi(x) = 0z>. 
Jm ' 
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Consequently Z' t is the exponential barycenter of i^(/xo) up to time T' n . This 
implies that T' n = T\ . 

On the set {T^ < r n } we replace time by time and time by 

T 2 n := r n A inf |t > T* , p{Z T , , Z t ) > ^ } , 

to prove the requested property up to time T% ■ With the same procedure we 
define a sequence of stopping times (T^)k>± converging almost surely to r n , 
such that Proposition 2.10 is true on [0, T™[. Consequently it is true on [0, T n [ 
as requested. 

□ 



2.5. The hyperbolic space example. Let M be a Riemannian manifold 
and let V be the Levi-Civita connection. Assume that M is convex for the 
Levi-Civita connection. Denote by p(x, y) the distance function between x 
and y which is abbreviated as p where there is no risk of confusion. For any 
two points x ^ y and 7 = j(x, y) the geodesic connecting them, a vector u £ 
T-rM has an orthogonal decomposition u = u L + u N where u L is in the direc- 
tion of 7(0). We may also use u L ( x ' y ^ and u N ( x ' y ^ when the points concerned 
need to be clarified. Denote by // x ,y ■ T X M — > T y M the parallel transport 
along the geodesic 7(2;, y). Set Ei(x,y)(s) = (// x ^^)Ei(x,y)(0) for some 
tangent vectors Ei(x,y)(0) on T X M such that (Ei(x,y)(0), . . . , E n (x,y)(0)) 

is an orthonormal frame of T X M . We shall fix Ei(x,y)(0) = ||^oj|| so that 

E 1 (x,y)(s)= T^fffv 

Since M is convex we can assume that the Ei(x, y)'s are chosen to be smooth 
in x and y in a neighborhood of some (xo>2/o) °^ diagonal. 

Denote by H m the hyperbolic space of dimension m and by G : HP™ x M m — > 
H m the smooth map which sends (x, y) to the center of the geodesic segment 
[x,y] -G{x,y) =exp a .(iexp" 1 y). 

Lemma 2.13. The map G:M m x HP™ -^HI m is harmonic. 

Proof. We need to prove that tr X7dG = 0. We do it with a symmetry 
argument. Observe that G(x,y) = G(y,x) for all x,y. We have 

(2.19) trVdG(x,y) =trVdG(y,x). 

On the other hand, the symmetry ip x ^ y in H m with center G(x,y) is an 
isometry which exchanges x and y. So 

tr VdG(y, x) = \xVdG{<p x , y {x),(p Xt y{y)) 

(2.20) 

= ( l f x ,y)*^^dG(x,y) = -trVdG(x,y). 



20 M. ARNAUDON AND X.-M. LI 

The identities (2.19) and (2.20) lead to tr X7dG{y,x) = 0. □ 

Proposition 2.14. For x ,y G H m set /j, = \{5 Xo + 5 yo ). Consider 
two independent Brownian motions Xt and Yt with initial points xq and 
yo, respectively. Set p t = p(X t ,Y t ) and define a measure fit by pt(A) = 
^(5x t {A) + ^Yt(^-))- Then Zt is a martingale in M. It satisfies the ltd equa- 
tion 

o*Z t = \{// Xu z t d v X t ) L + \uiy t ,z t d V Y t ) L 
2cosh{pt/2) 

where u = u L -\-u N with u L tangential (resp. u N orthogonal ) to the geodesic 
"f(Xt,Yt). Furthermore, we have locally 

1 1 m 

d v Z t = —E 1 (Z t ,X t ) dWl + E l (ZuX t ) dWl 

V2 V2cosh(p t /2) ^ 

/or Wj? independent real-valued Brownian motions. 

Proof. We have = G(Xt, Yt), where G is smooth and harmonic. Con- 
sequently, since (X, Y) is a Brownian motion in H m x M m , Zt is a martingale 
in EI m . We first look for a Stratonovich equation for Zt- All Jacobi fields 
along 7(3;, y) can be written as: J(s) = (a 1 + b 1 s)Ei(s) + Y^ l ™ =2 [ cosl[l {sp)a 1 + 
sinh(sp)6' i ]£'j(s), where Oj, hi G R and p = /9(x, y). In particular for u G T^H" 1 , 
v G TylH" 1 and J(u, v) the Jacobi fields along the geodesic j(x, y) with bound- 
ary values u and v, one has 

J(n,u)(s) = (u l {l -s) + sv^E^s) 

+ V f(cosh(sp) - sinh(sp) cothp)n* + smh ( gp ) V A EAs), 
~" 2 V sump J 

where u % = (u, £^(0))h, v l = (v, Ei(l))^ and (•, -)h denotes the scalar product 
in H = M m . Note that ^E^O) = u L ^ and v 1 E 1 (l) = v L ^' y \ This gives 

m 

j(«,0)(0) = -^(0) -Y / p(x,y)coth(p(x,y))u l E i (0) 

i=2 

(2.21) 

= _ u L(x,y) _ ^ y) coth ( p ( a . j j,))„^(*.v) ; 

j(o,« )(o) = ^iW+E . » «*^(o) 

^ smh(p(x,y)) 
swh(p(x,y)) 
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We have 

^j(SZ t ,0 Ftiy) )(0)dfi {y) 

= -i(^) i(z *' Xt) -i/9(Z t ,X t )coth(p(Z t ,X t ))(^) JV ( z *' x *) 
- \{5Z t ) L( - z ^ - ±p(Z t ,Y t )coth(p(Z t ,Y t ))(5Z t ) N ( z ^ 

= -(«5z t ) L(Zt ' Xt) -ip(x t ,y t )coth(i P (x t) y t ))( ( 5z t ) Ar ^ A ) ) 

since Xt, Yt and Z% lie on the same geodesic v L ^ Zt,Xt ^ = v L ^ Zt ' Yt > for any 
vector v. We shall abbreviate v L ( Zt ' Xt ^ =v L ( z t¥t) as V L Q n other hand 

j(0 Zt ,SF t (x))(0)dno(x) 
= ±j(0,5X t )(0) + ±j(0,5Y t )(0) 

^111 XV \ L i 1 P( Z tl X t) in XV \N 

= - 2 (l/x„zMt) + 2-j^^(//'.A«.) 

xv \ L i 1 gCfeijj) / // xv^ 
2 4sinh(±p(X t , Y t )J 

i xvm l i 1 P(Xt,Y t ) xv\ N 

+ 7jU/Y t ,Z t 0Yt) + - . — j v sA //Y t ,Z t SY t ) . 

2 4sinh(|p(At,y t )) 

Finally, applying Proposition 2.10(h) gives the Stratonovich equation 

<JZ t = \{//x t ,z t 8Xt) L + ^(//y t ,z i <5y) i 

+ \ cosh ( jp(X t) Y t ))((// Xt>Zt 6X t ) N + V/ Yt *6Ytf). 

Obtaining the Ito equation (2.14) from the Stratonovich equation is im- 
mediate since Zt is a martingale. 

For the local form of the Ito equation, by a localization procedure as in [8], 
Lemma 3.5, we only need to define W % locally. We do it with the formula 

(2.22) dWt = ±{Ei(Zt,X t ), // Xt ,z t d^X t + // YuZt d v Y t ) m , 

valid for (X S ,Y S ) in a small open subset of HP 71 x H m \ {diagonal}. The 
processes (W/) are martingales with quadratic variation t and therefore 
Brownian motions. The independence comes from the mutual orthogonality 
of EiiZuXtYa. □ 
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Remark 2.15. Note that if ujq is a point such that the geodesic (Xt(uJo)Yt(uJo)) 
converges, then Z t is close to the solution to dz t = -^E\(z t , X^^q)) dW^, 

which lives on the geodesic segment (X.(loo)Y.(uq)) oc if it starts there. In 
the next section we shall make precise this convergence. 

3. The barycenter of two independent Brownian particles in H. Con- 
sider the hyperbolic space EI = M. d of dimension d > 2 as an oriented man- 
ifold. Let dM be the visibility boundary in its visibility compactification, 
that is, the set of equivalence classes of geodesic rays. In the upper half 
space model dM is the union of the boundary hyperplane and the point at 
infinity. The visibility topology coincides with the usual topology when H is 
considered as a subset of M. d . 

Let X and Y be two independent Brownian motions in M with initial 
points Xq and Yq, respectively. Then they converge in EI with limits Xoo 
and in dM, respectively. We denote by (XY)oo C EI the random geodesic 
connecting the two limit points. By a random variable u on (XY)^ we mean 
u{ui) € (XFjoo^) for almost all u. 

Theorem 3.1. Let X% and Y-t be two independent Brownian particles 
in M and let Zt be their exponential barycenter. Denote by [i{Xq,Yq) the 
law of (Zt)t>o- Let u be an J-^ -measurable random variable on (XY)^ 
and let (T n ) n >i be a sequence of finite stopping times increasing to infin- 
ity such that Zt„ converges to u almost surely. Then for almost all ooq, 
fi(XT n (coo),Yr n (u>o)) converges as n goes to infinity to the law of a Brown- 
ian motion (zt)t>o of variance 1/2 on the geodesic (XY)^^^) with starting 
point u{u)q); that is, (z2t)t>o ^ s a Brownian motion on (Xy) 00 (o;o)- 

Proof. For the proof take the upper half space model 

M = {y=(y\y 2 ,...,y d )eR d ,y d >0}. 

The process (X t ,Z t ) is a diffusion as the image of the diffusion (X,Y) by 
the diffeomorphism (x, y) i-> (x, j(x, y){\)) on EI x EI. Set X t n = XT n +t, Y t n = 
Y Tn+t and Z£ = Z Tn+t . Then (X?, Z?) has initial condition (X Tn ,Z Tn ) which 
by assumption converges to (Xoo , u) . 

Since the hyperbolic Laplacian is given by Ae/(x) = ^(x d ) 2 A — ^ L x d d ( i, 
the hyperbolic Brownian motions can be written as the solution to the ltd 
stochastic differential equation 

dX t = X? dB t -\{d- 2)Xfe d dt, 

(3- 1 ) .1,1 w 

dY t = Y d dB' t -\{d- 2)Y t d e d dt, 

where (Bt,B' t ) is a Brownian motion in M? d and {e{}f =1 is the canonical 
basis in R d . 
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Set f(x l ,...,x d ) = sup(x d ,0). Then (X t ) can be considered to be the 
solution to 

dX t = f(X t ) dB t - \(d - 2)f(X t )e d dt 

on M. d . This equation has locally Lipschitz coefficients. We shall show that 
the process (Xt, Zt) also satisfies an equation which extends to M. d x EI with 
local Lipschitz coefficients. First note that since Xt and Yt are martingales, 

(3.1) yields 

d v X t = X d dB t , d w Y t = Y d dB' t . 
On the other hand, by Proposition 2.14 

d v Z, = \(Hx„z, d v X t ) L + \(I/y„z, <i v y,) 1 

+ a75) ((//x '' z ' <iVx,)A ' + (//y ' AdV> '' )N) - 

We have 

(3.2) dZ t = d w Z t - \Y(Z t )(d v Z t ,d v Z t )dt, 

where T(z) is the Christoffel symbol at z for the canonical connection in 
H. It is a smooth function, so we only need to prove that the equation for 
dy Zt extends as requested. For this plug d? Xt and d^Yt into the formula 
for d? Z t and observe that Y(t) = H(X(t),Z(t)) for H(x, z) = exp(2exp~ 1 z) 
to obtain 

d^Z t = \(llx u z t (^t) L{XuZt) 

+ // H (x u z t) ,zAmXt,Z t ) d dBl) L ^ z ^) 

// XuZt (X d dB t ) N ^ 
2cosh(p(Xt,Zt)) 

// H{ x t ,z t) ,z t (H(Xt, ZjfdB^fW^^ 
2cosh(ppQ,Z t )) 

Let w = (w 1 , . . . , w d ) be a vector in R d and let x n = (x\, . . . , x d ) be a se- 
quence in HI converging to x = (x 1 , . . . ,x d ~ 1 ,0) G dW \ {oo}. Set E(z,x) = 
p (l^ exp" 1 x, which extends smoothly to the set (Hxfl) \ {diagonal} con- 
sidered as a subset of i x R^. Furthermore, 

= lun a (-x d )(w,E(x n ,z)) m E(z,x n ) 

= lim (w,-(x d )~ 1 E(x n ,z))E(z,x n ) = -w d E(z,x) 
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since (x d ) 1 E(x n ,z) converges to e^, and similarly 

Moreover, for z £ EI, 

n ^°° 2cosh/}(a; n ,z) 

from the boundedness of the nominator and the convergence to +00 of 
cosh/o(x n , z). We conclude that each coefficient in the equation for (X t ,Z t ) 
smoothly extends over H x 1 and (Xt,Zt) is a solution to the system of 
equations of the following form: 

(3.3) dX t = f(Xf) dB t ~ ^f(X?)e d dt, 

(3.4) dZ t = a(X t , Z t ) d(B t , B[) + b(X t , Z t ) dt, 

where a(x, z) = a(x, z) and b(x, z) = b(x, z) if x d < and x = (x 1 , . . . , x d ~ l , 0). 
Since all coefficients are smooth on H x EI and (M. d \ H) x EI they are locally 
Lipschitz on x EL Consequently the system (3.3)-(3.4) has a unique so- 
lution and it does not explode for starting points in H x H as known. If, 
however, the starting point satisfies X d < 0, then X. = Xq and (3.4) reduces 
to 

(3.5) d v Z t = \E(Z t ,Xu)(-dB d + dB' t d ) 

whose solution starting from Z$ shall be denoted by (Zt(Xo, Zq)). Then 
Z2t{Xo, Zq) is a Brownian motion on the hyperbolic geodesic (XqZq). In 
particular, it does not explode. 

Consider the process (X n ,Z n ) as solution to (3.3)-(3.4) where (B t ,B' t ) 
is replaced by (BT n +t,B' Tn+t ), with starting point (XT n ,Zx n ) which by as- 
sumption converges almost surely to (X^ , u) . Since M. d x EI is diffeomorphic 
to R 2d we can apply Corollary 11.1.5 in [25] for M m -valued SDE and con- 
clude that the transition probabilities are Feller continuous and so by the 
Markov property the law of (X n ,Z n ) conditioned by (Xx n ,ZT n ) converges 
almost surely to the law of (A^, Z(X OQ ,u)). Hence the law of Z^t is that of 
a Brownian motion on the geodesic (Xoo,n), as requested. □ 

The rest of this section is devoted to the existence of sequences of stopping 
times (T n ) n>1 as in Theorem 3.1. 

First we note the following fact: If Xt and Yj are independent Brown- 
ian motions in EI, then for almost all u> the random geodesic {Xt{uj)Yt{u))) 
converges to {X(uj)Y {uj))oo uniformly on compact sets as t goes to infin- 
ity, that is, for any compact set K, lim^oo sup^ g(Xt y t)nAr p(z, (I OT ^oo)) = 0. 
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To see this we take the upper half space representation. We may assume 
K C {(x 1 , . . . ,x d ) G H, x d > e} for some e > 0. If z G n if, then z be- 

longs to a Euclidean tubular neighborhood of the half-circle (J7) w , with 
radius sup s>t max(||X s — -X^ll, — Since for z £ K the hyperbolic 

distance is smaller than the Euclidean distance divided by e, 

lim sup p(z, (XooYoo)) < lim - sup ( max(||X s - Xoo\\, \\Y S - 3^»||)) 
t_>00 a6(X t yt)nA- £ s>t 

= 0. 

Next for x,y,x G H with x ^ y, denote by (xy) the geodesic segment 
connecting them and by p(z,x,y) the orthogonal projection of z into the 
geodesic (xy). 

Lemma 3.2. Let (Xt), (Yt), (T n ) n>1 and u be as in Theorem 3.1. // 
(u n ) n>1 is a sequence of (pT n )- adapted random variables in H converging 
almost surely to an J 7 ^ -measurable random measurable u on (XY)oo, then 
]im n - too p(un,XT„,YT n ) = u almost surely. 

Proof. Since the projection to the convex set (XT n Yr n ) is 1-Lipschitz 
we have 

p(u,p(u n ,X Tn ,Y Tn )) 

< p(u, p(u, X Tn , Y T J) + p(p(u, X Tn ,Y Tn ), p(u n , X Tn ,Y T J) 

<p(u,(X Tn Y Tn )) + p(u,u n )^0 

following from the fact that for almost all to the geodesic (Xt(uj)Yt(u>)) con- 
verges to (I7)oo(w) uniformly on compact sets K of EI . □ 

Theorem 3.3. Let (Xt), (Yt), (Zt) and u be as in Theorem 3.1. There 
exists a sequence of finite stopping times (T n ) n >i increasing to infinity such 
that Zt„ converges to u almost surely. 

Proof. Let (u n ) n>1 be a sequence of (.T^) -adapted random variables 
in H converging almost surely to u. By Lemma 3.2 it is sufficient to prove 
that there exists a sequence of stopping times (T n ) n >i such that T n >n and 
Zt„ = p(u n , XT n ,Yr n ) almost surely. Conditioning with respect to T n , it is 
sufficient to prove that for every o G M and n > 1, there exists a stopping 
time T n >n such that Z Tn = p(o,X Tn ,Y Tri ). 

Let Pt = p(o,X t ,Y t ). Since Pt and Z t belong to the geodesic segment 
[Xt, It], Pt = Z t if and only if the signed distance D t = p(Pt,X t ) — p(X t ,Z t ) 
between Pt and Zt is zero. For the existence of the stopping times we only 
need to show that D t is recurrent. Set = p(o,X t ) and Rj = p(o,Y t ). By 
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the definition of Pt and the triangle inequalities, we have p(Xt,Pt) < R^ 
and p(P t ,Y t )>Bj -p(o,P t ). So 

D t = p(P t ,X t ) - \p{X t ,Y t ) = p{P t ,X t ) - \{p{X t ,P t ) + p(P u Y t )) 

= \p(X u P t ) - \p(P t ,Y t ) < \{R? - Rj) + \p{o,P t ) 

and similarly D t > \ {R^ — R[) — \p(o,Pt). Since limt->oo p(o, Pt) exists we 
only need to show R^ — Rj is recurrent. Note that 

dRf = dB? + coth R? dt, 

dRj = dBj + d -^- coth Rj dt, 

where (Bf ,E>J) is a planar Brownian motion. So for almost all ui, 

ini{R?{ui),Rj(uj))>\(d-l)t 
for sufficiently large time. Now 

rt 

R? - R( = B? - Bj + \ [d - 1) / (coth i?f - coth R Y s )ds. 

Jo 

But B^ — BY is recurrent and Jq°° (coth Rf — coth R\ ) ds exists since 

| coth^f - cothi^| = | coth^f - 1 - (cothi?J - 1)| 

2 2 

" e 2R* _ I ~ e 2Rj _ 1 

for large time s. Consequently R^ — Rj is recurrent and so is Dt- □ 

4. Barycenters of measures in a Cartan Hadamard manifold. Let M be 

a Cartan-Hadamard manifold with pinched negative curvatures k, —b 2 < 
k < —a 2 , for b > a > 0. Denote by M = M U dM its visibility compactifica- 
tion where dM is the set of equivalence classes of unit speed geodesies in M 
under the equivalence relation 

7i~72 limsupp(7i(i),7 2 (£)) < oo, 

t— >oo 

endowed with the sphere topology. See, for example, [2], page 22. Note that 
for each point z £ M and x E dM there is a unique unit speed geodesic 
in the equivalence class of x with initial point z, which we shall denote 
as {tp(z,x)(t), < t < oo}. For z,y £ M, z ^ y, we shall also denote by 
{(p(z, y)(t),t > 0} the unit speed geodesic satisfying ip{z, y)(0) = z and cp(z, y)(p(z 
y. In other words, 

^P(z,y)(t) =~f(z,y)(t/p(z,y)) \/z,ye M satisfying z^y. 
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Let 7r : TM — > M be the canonical projection and let 

: SM = {v£ TM, \\v\\ = 1} — ► dM 

be the map which sends 6 to [exp 7r( - 6 ,)(t^)] t >o S <9M. The map $ = (tt,4>) 
from SM to M x <9M is a homeomorphism. In fact, = ip(z,x)(0) 

(see, e.g., [6], Propositions 2.13 and 2.14). 

For o E M denote by (ip ,x, % G <9M) the family of Busemann functions: 

(4.1) Vo,x(y) = lim [p(y,¥>(o,x)(*)) -t]. 

t — >oo 

These functions are characterized by i/j ,x(o) = and grad^ Vo.z = — 3? u/j a?) = 
—<p(y, x)(0), any y S M. We write V'x for ip ,x if there is no risk of confusion. 
Denote by A4 1 (X) the set of probability measures on a topological space X 
endowed with the Borel cr-field. Set 

(4.2) ^{z) = ^ 0jtl (z)= [ ^ 0jX (z)dfi(x), fieMi(dM). 

JdM 

Then 

(4.3) grad z W = - / ip(z,x)(0) d(i(x). 

JdM 

A solution to gradj.?/^ = is called a Busemann barycenter of \i. 

The aim of this section is to show that for discrete probability measures 
of finite support transported by a suitable random flow, the exponential 
barycenter fa converges to the Busemann barycenter of the limiting measure 
on the boundary. 

Lemma 4.1. Let fx be a probability measure on M with corresponding 
vector field M — > TM defined by (2.1). We have 

(VuH^u) < -ae(z,fi)\\u\\ 2 , u£T z U, 

where 

(4.4) e(z,/j,)= min / p(z, x) sin 2 (w, <p(z, x)(0)) d/j,(x). 

w£S z M Jm 

Proof. We follow the notation of Section 2.3 and note that 

W u H il = f j(u,O x )(0)dn(x), 
Jm 

where X £ T X M is the zero vector. Write u = u L + u N where u L = u L<yZ,x ^ is 
colinear to expj 1 x and u N = u N ( z ' x ' is its orthogonal complement. Then 

j(u,O x )(0) = j(u L , 0,5(0) + J(u N , X )(0) = -u L + j(u N , X )(0), 
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sum of two orthogonal vectors. Write J(t) = J(u N ,0 x )(t) and f(t) = \\J(t)\\. 
Then f{t) = jjftj\\(J(t)J(t)). In particular, f(0) = ■Aq(u N i J(0)) and fur- 
thermore 

f"(t) = r 3 (t)(|| J(t)|| 2 || J(t)|| 2 - (J(t)J(t)f + \\J(t)\\ 2 (J(t),J(t))) 

>r\t)(j(t),j(t)) 

= -r 1 (t)(J(t),R(J(t),j(t))j(t))>a 2 p\z,x)f(t), 

where "y(t) = "y(z, x)(t) for t 6 [0, 1] and a is the upper bound of the sectional 
curvature. Note that /(0) = \\u || and /(l) = 0. By comparison with the 
solution to g"(t) = a 2 p 2 (z,x)g(t) with the same boundary conditions, we see 

/'(0) <g'(0) = -ap(z,x)coth(ap(z,x))\\u N \\ < -ap(z, x)\\u N \\ 

which leads to (u N ,j(0)) < -ap(z, x)\\u N \\ 2 . Finally 

{V u H IM u)= j (-\\u L \\ 2 + (j(0),u N ))dp(x) 
Jm 

<l (-h L \\ 2 -ap(z,x)\\u N \\ 2 )dp(x)<- f ap(z,x)\\u N \\ 2 dp(x) 
Jm Jm 

= -a\\u\\ 2 l p(z,x) sin 2 K**) <-IMN^)- D 

Jm lj 

Lemma 4.2. Let p be a probability measure on dM . Then 

(4.5) VdTp fl (u,u)>aa(z,p)\\u\\ 2 Vu£T z M, z£M, 
where 

(4.6) a(z,p)= min / sin 2 (w, <p(z, x)(0)) dpix). 

w£S z M Jqm 

Proof. Let z £ M and x G <9M. Set 

Vv(-) = p(-> ^(^o, z)(*)) - t. 

By Proposition 3.1 in [13], as t goes to infinity, ('ipt,x, gradf/^z, V grad^x) 
converges to (V>z , grad -i/^ , Vgrad'i/'x) uniformly on compact sets. In fact, the 
proof of the same proposition shows that the convergence is uniform in x. 
In particular, if we set 

AA Z ) = / i>t,x(z)dp(x), 

JdM 

then 

grad z W = - / ^>(z,ip(z ,x)(t))(0)dp(x), 

JdM 

V u gradVv = - / V u (p(-,<f(z ,x)(t))(0)dp(x), 

JdM 
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(i/^grad^, Vgrady>£ iA1 ) converges to (^ M ,grad^, Vgrad^ M ) uniformly 
on compact sets and the convergence is uniform in /i. 

Let ^ = ^(z ,-)(i)*(/u). Since j(x,y)(0) = p(x,y)(o)tp(x,y)(0), 

H IH {z)= I j(z,y)(0)dp t (y) 

JM 

p(z, 99(20, x)(t))tp(z, ip{zo,x)(t))(0) dp{x). 

dM 

If u G T ZQ M, then d u p{-, <p(zo, x)(t)) = —(ip(zo,x)(0),u) and 

V u i^ t = / p(^b,<^(^,a;)(t))V, t ^(-,^(2o,x)(t))(0)d/i(a:) 

+ / d u p(-,ip(z Q ,x)(t))(p(z ,ip(z ,x)(j;))(0)dp(x) 
JdM 

tV u (p(-,(p(z ,x)(t))(0) dp{x) 

dM 

(<p(zo, x)(0),u)ip(zo,x)(0)dp(x) 

dM 



= -tV u gradip t ,ti ~ / (<p(z Q ,x)(0),u)ip(z Q ,x)(0)dp(x) 

JdM 

using p(zo, ip(zo, x)(t)) = t and for all s > 0, (p(zo,<p(zo, x)(t))(s) = ip(zo,x)(s) 
Consequently 

(V u Hu t ,u) = -f(V u grad^ „,*u) - / (^(z ,2;)(0),u) 2 (i / u(2;). 

This together with Lemma 4.1 gives 
(V u grad ipt,n,u) 

> ^(zo,Pt)\\u\\ 2 - - [ {(p(z Q ,x)(0),u) 2 dp(x) 



t t JdM 

d f 

>tII u I| 2 min / p(z ,y)sm 2 (w,<j>(z Q ,y)(0))dp t (y) ~ -\\u\\ 2 
t w£Sz MJm t 

>7lM| 2 m ir \ f / P(*b, a:) (*)) 

t w£S ZQ M Jqm 

x sin 2 (w,ip(z , tp(z ,x)(t))(0))dp(x) - ^\H\ 2 

= a||u|| 2 min / sin 2 (w, ip(zo, x)(0)) dp(x) ||u|| 2 

w&Sz MJqm t 

= aa(zo,/u)|M| 2 — yll u l| 2 - 
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Taking t to infinity and using the convergence of V gradt/^ to Vgrad-0^, 
we obtain V 'di^^{u,u) > aa(zo, //)||u|| 2 , as desired. □ 

Let IA be the subset of M\{dM) containing discrete measures with no 
atoms of weight greater than or equal to 1/2. Clearly the function a(z,fi) is 
continuous onMxW, using the weak convergence topology on U. Moreover, 
it is strictly positive on M x U: for every z G M, the set {ip(z,x)(0), x G 
supp(/x)} contains at least three different vectors, and this implies that no 
w G S Z M is colinear to all of them. The positivity then follows from the 
compactness of S Z M. 

Proposition 4.3. For fi£U there exists a unique z G M such that 
(4.7) grades- / ' (f(z,x)(0) d/i(x) = 0. 



IdM 

Denote the solution by G{jj). Then the map G:IA — > M is continuous. 

Proof. The existence and uniqueness are well known in the case that \i 
is a continuous measure. In fact the uniqueness follows from Lemma 4.2 since 
a(z, fi) > and ipn is strictly convex. For the existence we only need to show 
that there is a geodesic ball B{o,T) C M of radius T > on which grad^ 
points outward the boundary and therefore has a zero inside the geodesic 
ball. For T > 0, take y £ dB(o, T). Let 7 G dM be the point corresponding to 
the geodesic ray 7(0, y) and let B e (fy) C dM be the set of points whose angle 
with 7(0, y) is smaller than e using the sphere topology. Choose £0 > so 
that Co = SM]? x ^dM M-Seo( x )) < 1/2 which is possible due to the compactness 
of the sphere at infinity. Let x G dM, e the angle between (p(0, y)(0) and 
ip(0,x){0), and a the angle between ip(y,x)(0) and (p(y,o)(0). By comparison 
with a manifold with constant curvature —a 2 , we have 

cos e cos a + 1 2 
sin a < . ; < 



sinecosh(aT) sinecosh(aT) 
and 

<0(y,x)(O),0(y,o)(O))>l 



sin e cosh (aT) 
Consequently 

-(grad y ^7i,<p(y,o)(0)) 



+ / ){<p(y,x)(p),<p(y,o)(P))Mx) 

— f fi - . 2 \ 2 , T + (-1) • KB £0 m 

JdM\B eQ (*i) \ sure cosh [al)J 
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= f 1 " ■ 2 \ 2( T . ) (1 " M(£ £0 (7))) - mOMxO 
\ sin e cosh (ai J / 

4 

> 1 k 27 7 ~~ ^co 

sin e cosh (aT) 

>0 

when cosh 2 (aT) > 4/ sin 2 e(l — 2co). Since (fi(y, o)(0) is normal to the bound- 
ary dB(o,T) and pointing inwards, this proves that grad^^ is pointing 
outwards. Consequently there exists z G B(o,T) such that grad., Vv = 0, so 
G is well defined. 

For the continuity of G note that grad y W ^ s continuous with respect to 
(y,/i) since grad y ipt,ti, which is continuous in (y,t,fi), converges to grad^V 
uniformly on y in compact sets and uniformly in fj,. Let (/i n )n>i be a se- 
quence of elements of U converging to /i G U. Set z n = G(fi n ) and z = G(fi). 
From the convergence of the sequence (fi n )n>i to /i, we can choose the 
same £o, cq and T for all \x n and /i. Consequently, z and all the z n be- 
long to B(o,T). Furthermore, the function a in Lemma 4.2 is continuous, 
so we have ao — inf{a(y,/x n ) : (y,n) G B(o,T) x N} > and consequently 
( V u grad , jUn , u) > a«o ||ii|| 2 for all y G B(o,T) and u G T y M. Let z n (t)t>o be 
C 1 paths in M with z n (0) = z and z n (t) = — gradvp^ n (z n (t)) . Then by dif- 
ferentiating || grad 4> ti n {z n {t)) || 2 with respect to t we see || gradip IJiri (z n (t))\\ < 
e _aao *|| grad^^^O))!!- It follows that z n = z n (oo) = lim^oo z n (t) and the 
length of the curve from z = z n (0) to z n is smaller than || grad-i/Vn ( z )\\/( aa o)- 
Since gr&dip^ n (z) — > grad'0 M (z) = 0, linin^oo p(z, z n ) = 0. So G is continuous. 
□ 



Example 4.4. Let EI C C be the Poincare upper half plane. The bound- 
ary of EI is the real line M, plus one point at infinity. For n > 2 and x\ < 
■ ■ ■ < x n in E set = ^ Y^j=i &xj ■ If n = 3, the angles between the vectors 
ip(G(^),Xj)(0) are ±2-7r/3. Considering first the case x 3 = oo and then the 
general case via homographic transformations, one finds 

_ xix 3 + x 2 x 3 - 2x ± x 2 + iV3x 3 (x 2 - xi) 
2x 3 - x\ - x 2 + iy/3{x 2 - x\) 

For n = 4, the situation is even simpler: G(/i) is the intersection between the 
hyperbolic geodesies {x\x 3 ) and (X2X4). 

Next we consider a Brownian flow Ft in M, that is, a semimartingale flow 
with characteristic (a(t, x, y,Lo),0, t) as defined in Proposition 2.10, such that 
for every x G M, Ft(x) is a Brownian motion [which is equivalent to saying 
that a(t, x,x,uj) = J2i=i e i ® where (ej)i<.;<d is an orthonormal basis of 
T X M\. We furthermore assume that Ft is unstable, that is, whenever y\ and 
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U2 are two distinct points in M, the distance between Ft{y\) and ^(2/2) 
converges to +00 as t goes to infinity. Typical examples of unstable flows 
are given by isotropic Brownian flows in the hyperbolic plane with positive 
first Lyapounov exponent (see, e.g., [24]). The following result is immediate. 

Proposition 4.5. Let F t be an unstable Brownian flow in M . Take 
e > 0, y\^V2 and let X\ = F t {yi), i = 1,2. Then for almost all to there is 
t(u>) such that p{X} ,Xf) > ((d — l)a — e)t whenever t > t(u). 

Proof. For x\ 7^ 22 and u = (ui,u%) G Tr Xl X2 -\(M x M) write Ui = Vi + 
Wi, the orthogonal decomposition with Vi tangential to the geodesic (x±X2)- 
We have by the Hessian comparison theorem, Vdp(u, u) > atanh(^)(||u>o|| 2 + 
H^i || 2 ) (see, e.g., [23] Lemma 1.1.1, where a similar calculation is done in 
positive curvature). As a consequence, the drift of T t := p(X^,Xf) is larger 
than (d — l)a tanh(^y i ) dt. When t is sufficiently large, then is large by 
instability of the flow, so its drift is close to (d— l)adt. On the other hand, 
since p is 1-Lipschitz, the process has a bracket satisfying d(T, T)t < 2dt. 
So we can conclude that for e > and t large, Tf > {{d — l)a — e)t. □ 

Proposition 4.6. Let Xt and Y% be two M -valued continuous func- 
tions converging to X^ and Y^ on dM, respectively, as i — > 00. Suppose 
furthermore that lim t _ >00 (p(o, X t ) — £t)/t = 0, lim^ 00 (/?(o, Y t ) — £t)/t = 
and p(Xf,Yt) > a't for t large, for some constants a',£ > and o £ M . Then 
Xoo 7^ ^00 • 

Proof. We only need to establish that if Pt is the orthogonal projection 
of Yf on the half-geodesic [oX t ) := {y?(o, X t )(s), s > 0}, then lim^oo p(Pt, Y t ) = 
00. We will prove that, for t large, p(Pt,Y t ) > ^^-t. When P t = o, this is 

clear. Let us consider the case where Pt 7^ o. Let t be a time such that 
p(o,P t )>p(o,X t ). Then 

p(o,Y t ) > p(o,P t ) = p(o,X t ) + p(X t ,P t ) > p(o,X t ) + p(X t ,Y t )-p(P t ,Y t ), 

which implies 

p(Pt,Y t ) > p(o,X t ) - p(o,Y t ) + p(X t ,Y t ). 

This clearly implies p(Pt,Yt) > when t is large. Now let t be a time such 
that p(o,P t ) < p(o,X t ). We have 

p(P t ,Y t )>p(o,Y t )-p(o,P t ), 

p(Pt,Y t )>p(X t ,Y t )-p(P t ,X t ). 

Adding the two together gives 

2p(P t , Y t ) > p(o, Y t ) + p(X t , Y t ) - p(o, X t ) ~ p(X t , Y t ) > a't 



BARYCENTERS OF MEASURES 



33 



which again proves that p(Pt,Yt) > when t is large. Thus ^ Yx. □ 

A stochastic process Xt is said to satisfy the law of large numbers with 
limit £ if limj^oo jp(o, Xt) = t. It is known that all Brownian motions satisfy 
the law of large numbers with nonzero limit if (a) M = M. d , (b) M is the 
universal cover of a compact Riemannian manifold with negative curvature 
(see, e.g., [19]), or more generally, (c) M is a Cartan-Hadamard manifold 
such that for some point o G M, Ap(o, — ) converges to a negative constant 
as r goes to infinity. 

Note that if p(X l (t),y) are more or less all of the same length independent 
of i, that is, p(X i (t),o) = f(t) + R^t) with R}(t) small compared to f(t), 
then the minimizer of YaL\ p{X l {t),y) is close to (or the same as) that of 
Y^aL 1 P 2 (X l (t), y). This is the consideration behind the following theorem. 

Theorem 4.7. Let p be a discrete probability measure on M with finite 
support and no weight greater than or equal to 1/2. Suppose Ft is an unstable 
Brownian flow satisfying the conditions of Proposition 4.5 and such that 
lim^oo 7/9(0, Xt) = £ where £ > 0. Denote by Zt the exponential barycenter 
of the pushed forward measure pt = Ft{p)- Let p^ be the measure F ao (p) on 
dM . Then limt— ,00 Zt = G(poo) almost surely. 

Proof. The measure p^ is carried by a finite set, and by Propositions 
4.5 and 4.6, for almost all uj, F oc (yi,uj) 7^ F oc (y2,uj) if y\ 7^ y2- This implies 
Poo £ U. For pE M fixed denote by @ p (q) the point on dM determined by 
the geodesic (f(p,q), that is, Q p (q) = 4>((p(p,q)(0)). Then F t {x) induces a 
measure p p> t on dM by 6 P (-) and p Pt t — > ^00 • Furthermore, by continuity 
of G, Proposition 4.3, lim^oo G(p Ptt ) = G(p OQ ). 

Define 

-p(o,F s (y))-£ 



Re 



dp(y). 



1 11 

Since p has finite support, lim,,-,^ R s = almost surely. 

Set Z[ = G{pt)- By definition fQ M (p(Z' t ,x)(0)dpZ o ,t(x) = 0, so 

H tls {Z' s )= [ j(Z' s ,F s (y))(0)dp(y)-£s f <p(Z' s , x)(0) dp , s (x) 

JM J8M 

= f (j(Z' s ,F s (y)M-^(Z' s ,F s (y))(0))dp(y) 
JM 

+ £s { (ip(Z' s ,F s (y)m-^Z' sl Q o (F s (y)))(0))dp(y) 



IM 

(p(Z' s ,F s (y)) - £s) v(Z' s ,F s (y))(0) dp(y) 

M 

+ £sf (^(Z' s ,F s (y)m ~ ^Z' s ,Q o (F s (y)))(0))dp(y). 

JM 
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Set R' S = £J M \\<p(Z> s ,F s (y))(0) - ip(Z' s ,@ (F s (y)W)\\ Mv)- By the conver- 
gence of Z' s , the angular convergence of Brownian motions and the con- 
tinuity off diagonal of the map (z,x) i— ► ip(z,x)(0) from M x M to SM, 
lim^oo R' s = almost surely. We have 

WHu. (Z' s )\\ < p(o, Z' s ) + R s s + sR' s = p(o, Z' s ) + s(R s + R' s ). 

On the other hand, for y £ M, let R s (y) = \ ^p(o,F s (y)) - 1\. Note that for 
z £ B(Z' S , 1), 

p(z,F s (y)) > p(o,F s (y)) ~ P (z,o) > is - R s (y)s - [p(z, Z' s ) + p(Z' s , o)] 
>£s-R s (y)s-l-p(Z' s ,o). 
Consequently, by Lemma 4.1, if u 6 T Z M and z £ B(Z' S , 1), 
(V u H„ s ,u) 

<-a||u|| 2 mm / p(z, sin 2 (w, (p(z, F s (y))(0)) dp(y) 

w£S z M J m 

<— a\\u\\ 2 \£s min / sin 2 (w, ip(z, F s (y))(0)) dp(y) 
\ weS z M Jm 

- R s s- p(o,Z' s ) - 1 

= -a\\u\\ 2 (£sa(z,p Z)S ) - R s s - p{o,Z' s ) - 1), 

where a is the continuous and positive function defined in Lemma 4.2. 
Since p ZjS converges to F oc (p) and Z' s converges in M to G(/Uoo), there 
exist C\{uj) > and s(u>) > such that for s > s(lo), 

(4.8) {V u H^ s ,u) < -CJufs VzeB(Z' s ,l), ueT z M. 

Next let (3 s (t), t > 0, be the C 1 path starting from Z' s and such that ^(5 s {t) = 
—H^^fis^t)). Then with an argument similar to the end of the proof of 
Proposition 4.3, (4.8) implies that for every t > smaller than the exit time 
from B(Z'„,1), 

\\H, a (Ps(t))\\<e- Cls \\H, s (Z' s )\\ 

and 

s 

(4 ' 9) p(o,^) + S (fl s + #J = p(o,Z^)/s + + 

Cis Ci 

Since the right-hand side converges to as s goes to infinity, we see that for 
s large, (3 s {t) stays in B(Z' S , 1) for all t > 0. Moreover, f3 s (t) converges as t 
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goes to infinity to j3 s (oo) = Z s which is the only point in the manifold where 
H Us vanishes. Prom (4.9), we get 

p(Z s , Z s ) < — . 

The right-hand side goes to 0, so Z s converges to G(/Xoo). □ 



Remark 4.8. When M is the hyperbolic space i m and a has second- 
order space derivatives almost surely bounded, Theorem 4.7 generalizes to 
a discrete measure ji with compact but not necessarily finite support. The 
proof is the same; the only difference is that one has to find another argument 
for the almost sure convergence of R s to 0. The new argument is as follows. 
In an exponential chart ^ based at o, let (a', b' , t) be the characteristic of the 
flow. Then a' has second-order space derivatives almost surely bounded and 
b' = |A^/ has first-order space derivatives almost surely bounded. Since the 
chart is centered at o, the distance to the origin is the same in the chart and 
in the manifold. Using [5], Theorem 2.1, we can say that almost surely, for all 
y G supp(/i), R s (y) = \ ^p{o, F s (y)) — £\ is bounded by a constant depending 
only on the support of fj,. Since R s (y) converges almost surely to as s goes 
to infinity, by dominated convergence, R s goes to 0. 
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